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Abstract 

Pulsed  power  hydrodynamics  experiments 
typically  involve  multi-megajoule  capacitor 
banks,  imploding  thin-walled  cylinders  of  solid- 
density  metal  by  the  electromagnetic  force 
associated  with  axial  currents  in  excess  of  ten 
megamperes.  Much  of  this  activity  operates 
with  technology  developed  for  other 
applications  (e.g.,  soft  X-ray  generation),  which 
can  be  significantly  larger  than  actually  required 
for  the  scientific  objectives  of  such  experiments. 
This  is  particularly  true,  if  imaging  diagnostics 
(e.g.,  proton  radiography)  with  greater  spatial 
and  temporal  resolution  become  available, 
pennitting  implosion  experiments  of  smaller 
scale  size.  The  present  paper  develops  scaling 
relationships  for  electromagnetically-driven 
implosions  and  indicates  the  substantial 
reduction  in  system  energy  with  reduced  scale 
size. 

I.  INTRODUCTION 

Pulsed  power  techniques  provide  excellent 
opportunities  for  investigating  matter  under 
extreme  conditions.  One  category  of  technique 
uses  multi-megampere  currents  to  implode 
solid-density,  hollow  cylinders  ( liners )  in  the 
manner  of  a  z-pinch  for  which  axial  current 
interacts  with  its  own  magnetic  field  to  create  a 
radial  implosion  force.  Such  implosions  readily 
provide  speeds  of  several  km/s,  and  are  capable 
of  creating  strong  shocks  and  damage  in 
converging  geometries  that  provide  stresses  in 
both  radial  and  azimuthal  directions. 
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The  precision  of  electromagnetic  implosion  can 
be  better  than  1%  in  space  and  time,  pennitting 
careful  diagnosis  of  physical  phenomena  and 
close  comparison  with  computer  calculations. 

A  principal  difficulty  in  recent  years  has  been 
supporting  facilities  to  provide  the  necessary 
multi-megampere  currents  for  liner  implosion. 
In  addition  to  the  cost  and  size  of  the  pulsed 
power  facility,  the  consequences  of  the  liner 
implosion  itself  favor  lower  energy  events  in 
terms  of  reduced  collateral  damage,  smaller 
replacement  parts  and  easier  connection  to 
diagnostic  systems.  Recent  advances  in  proton 
radiography  as  a  high  resolution  diagnostic 
technique  and  continued  interest  in  experiments 
with  environmentally- challenging  materials 
separately  call  for  the  development  of  pulsed 
power  and  liner  implosion  systems  that  can  be 
brought  to  special  locations  (e.g.,  proton  beam- 
line).  We  explore  here  the  necessary  scaling 
relations  for  reducing  the  system  size  while 
retaining  important  intensive  values,  e.g.,  liner 
speed,  compression  ratio  and  energy  density. 

II.  BASIC  FORMULATION 

To  provide  true  scaling  between  pulsed 
power-driven  liner  implosions  as  we  reduce  the 
system  size,  we  should  attempt  to  achieve  the 
same  implosion  trajectory,  speed  and 
thermodynamic  state.  It  is  therefore  useful  to 
normalize  the  equations  for  the  liner  and  driver- 
circuit  behavior  in  terms  of  characteristic  values 
and  dimensionless  parameters.  We  start  with  a 
simple  LRC-circuit  driving  an  imploding-liner 
load  of  increasing  inductance  and  resistance  as 
depicted  in  Fig.  1. 
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Figure  1.  Schematic  of  simple  LRC-circuit  driving  liner. 

The  equations  for  this  system  in  dimensional 
fonn  are: 

[Lt  +  Lp  +  L(t)]  dJ/dt  =  V  -  J(dL/dt  +  R)  (la) 


dV/dt  =  -  J/C  (lb) 

dn/dt  =  (J2/2m)dL/dr  (lc) 

dr/dt  =  -  u  (Id) 

dL/dt  =  u  dL/dr  (le) 

dQ/dt  =  r|0(J/A)2  [  1  +  PQ  ]  (If) 

R  =  R0[1  +  PQ]  (lg) 


With  reasonable  transparency,  we  denote  the 
liner  mass  by  m,  inductances  by  L,  resistance  by 
R  and  capacitance  by  C.  Thus,  LT  and  Lp  are  the 
inductances  of  the  transmission  plates  and 
parasitic  inductance  of  the  capacitor  and  switch, 
respectively;  the  total  initial  inductance  is  L0  = 
Lt  +  Lp.  The  resistance  displayed  in  Eqn.  lg  is 
only  that  of  the  liner  itself,  but  will  be 
augmented  later  by  additional  (constant)  series 
resistance  of  the  plates,  switches  and 
connections.  The  resistivity,  rj,  and  resistance, 
R,  of  the  liner  increase  linearly  (from  initial 
values  denoted  by  the  subscript  ‘o’)  with  the 
energy  per  unit  volume,  Q,  deposited  by 
resistive  heating.  This  is  adequate  for  the  liner 
remaining  in  the  solid-state,  which  is  consistent 
with  our  neglect  of  many  complex  processes 
captured  more  appropriately  in  a  fully-posed 
MHD  computer  code.  In  deference  to  those 
who  use  ‘I’  for  impulse,  we  denote  the  total 
current  by  J. 


To  focus  quickly  on  issues  of  liner  size,  we 
have  simplified  the  formulation  to  a  ‘thin-liner 
approximation’  in  which  the  current  density  is 
uniform  in  the  liner  thickness.  The  thickness 
increases  during  implosion  to  conserve  the 
cross-sectional  area,  A,  in  the  manner  of  an 
incompressible  fluid,  but  we  retain  a  unifonn 
current-density  approximation  nevertheless. 
This  approximation  permits  a  first-order 
estimate  of  the  liner  heating  in  the  solid-state 
and  may  be  checked  a  posteriori  to  verify  that 
we  have  not  melted  the  liner  completely. 
Indeed,  we  should  expect  that  high  magnetic 
field  values  on  the  outer  surface  of  the  liner  will 
result  in  melting  there.  Prior  to  that,  however, 
we  would  have  more  than  doubled  the  electrical 
resistivity  of  the  material,  which  corresponds  to 
nonlinear  magnetic  diffusion,  tending  to  justify 
the  uniform  current-density  approximation  and 
spreading  the  heat  load  through  the  liner 
thickness. 

Nonnalization  of  Eqns.  1  follows  by 
introduction  of  characteristic  values  for  all  the 
variables.  In  some  cases,  the  choice  of  a 
characteristic  value  is  evident  a  priori,  such  as 
use  of  the  initial  value  of  voltage,  V0,  on  the 
capacitor  to  scale  subsequent  voltages,  and  the 
initial  liner  radius,  r0,  to  measure  the  radius 
further  in.  Reference  to  a  simple  LC-circuit 
suggests  also  that  the  characteristic  values  for 
time  and  impedance  are  (L0C)  and  (Lq/C)  ", 
respectively.  These  choices  imply  a 
characteristic  values  for  current  of  V J(LJC)xa 
and  for  inductance,  L0.  We  then  have 
dimensionless  variables: 

x  =  t/(L0C)1/2 ,  v  =  V/V0  ,  j  =  J/lVo/tWQ1'2] 

k  =  L(t)/L0  ,  0  =  R0/(L0/C)1/2  ,  q  =  [1  +  pQ  ]  (2) 

co  =  u/uc  ,  p  =  r/r0 

where  uc  is  a  characteristic  speed  not  yet 
specified.  Substitution  of  these  variables  into 
Eqns.  1  provides  dimensionless  differential 
equations: 

dj/dx  =  [  V  -  j  (  co/ p  +  0C  +  0q)  ]  /  (1  +  k )  (3a) 
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dv/dx  =  -j  ,  dk/dx  =  co/p 

(3b,  c) 

dco/dx  =  rj2/p  ,  dp/dx  =  -  Aco 

(3d,e) 

dq/dx  =  Kqj2 

(30 

In  this  process,  we  find  that  a  useful  value  of 
characteristic  speed  is: 

uc  =  Ldo/(L0/C)1/2 

(4) 

where  Ld0  is  the  initial  value  of  the 
inductance  gradient  dL/dr.  We  also  define  five 
dimensionless  parameters: 

Dynamic  parameter 
r=[CV02/2m]/[(L0/C)1/2/Ld0]2 

(5a) 

Inductance  parameter 

A  —  L0/Ld0i*o 

(5b) 

Liner  and  circuit  resistance  parameters 

0  =  Ro/(Lo/C)1/2  ,  0c  =  Rc/(Lo/C)1/2 

(5c, d) 

Liner  heating  parameter 

K  =  pn0CV02/[A2(L0/C)1/2] 

(5e) 

while  preserving  the  normalized  liner  and  circuit 
behaviors: 

Initial  thickness 

§0  =  (n0/|i)  [  2(3/7  r  /  K0  ] 1/2  (6a) 

For  aluminum,  r|0=  2.7x1 0^12-m.  p=2.1xl0"9m3/J,  and  p  = 
2.7xl03kg/m3,  so  with  T=0.75,  0=0.001,  and  K=2,  80  = 
1.4  mm.  Note  that  this  is  independent  of  scale  size. 

Inductance 

L0  =  (ph/2n)  A  (6b) 

=  20  nH  ,  for  A  =  2  and  h  =  5  cm 

Capacitance 

C  =  2pp0rAro2/(|ih/27t)K  (6c) 

=  2126  |if ,  for  r0=  5  cm,  h  =  5cm 

Characteristic  impedance 

Zc  =  (L0/C)1/2  (6d) 

=  3.1  m£2 

Characteristic  time 

tc  =  (L0C)1/2  (6e) 

=  6.5  jus , 

for  x  =  1.6  (rQ/r  =  2),  t  =  10.4  (is 

Liner  speed  at  above  time,  with  ta  =  0.4 

u  =  co  Z  J  Ldo  (6f) 

=  6.2  km/s 


Specification  of  these  dimensionless  parameters 
provides  nonnalized  solutions  for  the  circuit  and 
liner  behavior  as  shown  in  Figures  2  and  3. 
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Figure  2.  Normalized  solution  for  current,  voltage  and 
liner  radius,  T  =  0.75,  A  =  2,  0  =  0.001,  0C  =  0.2,  K  =  2. 


Initial  Stored  Energy 

W0  =  2jir0h  (rjo/pit  K p  T/2p  ]1/2/03/2  (6g) 

=  10.5  MJ  ,  corresponding  to  V0  =  99  kV 
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Figure  3.  Normalized  solution  for  resistance  increase  and 
liner  speed,  Y  =  0.75,  A  =  2,  0  =  0.001,  0C  =  0.2,  K  =  2. 


III.  SCALING  RELATIONS 

For  any  set  of  the  dimensionless  parameters 
(Eqns.  5),  we  can  return  to  dimensional  units  in 
terms  of  the  initial  radius  and  length  of  the  liner, 


If  we  change  the  liner  dimensions  to  r0  =  2  cm 
and  h  =  1cm,  the  various  dimensional  quantities 
change  (except  for  the  liner  thickness).  Of 
particular  concern  is  the  decrease  of  initial 
circuit  inductance  to  4  nH,  along  with  a 
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reduction  of  stored  energy  to  0.85  MJ. 
Competition  with  minimum  values  for  parasitic 
inductance  of  the  capacitors  and  switches  thus 
begins  to  intrude  on  the  system  design. 
Furthennore,  the  peak  magnetic  field  on  the 
outer  surface  of  the  liner  increases  from  about 
one  megagauss  for  r0  =  5  cm  to  1.65  MG  for  the 
smaller  liner.  In  the  linear  diffusion  regime,  the 
surface  temperature  scales  as  the  square  of  the 
magnetic  field  increasing  the  challenge  of 
controlling  the  liner  material.  Both  the  magnetic 
field  and  the  electric  field  scale  as  r0"  in  this 
fully-diffused  model,  (and  are  independent  of 
liner  length). 

IV.  OTHER  CIRCUITRY 

Transformer  techniques  can  mitigate  the 
effects  of  decreasing  circuit  inductance  with 
reduced  size  by  allowing  a  better  match  between 
the  primary  storage  circuit  and  the  relatively 
low  inductance  change  associated  with  liner 
implosion  (~  several  nH).  A  system  in  which 
each  turn  of  the  transfonner  primary  has  an 
opening-switch  has  the  best  electrical 
perfonnance,  but  can  be  cumbersome.  A  simple 
current-step-up  transformer  coupled  directly  to 
the  liner  load  (Fig.  4)  may  be  adequate. 


The  dimensionless  parameters  are  similar  to 
the  direct-drive  case,  but  referenced 
appropriately  to  the  primary  or  secondary 
inductances.  We  also  need  additional 
parameters,  such  as  the  usual  coupling 
coefficient,  k  =  M/(LPLs)1/2,  and  nonnalized 
uncoupled  inductances 

k,  =  Lj/Lp,  k0  =  L0/Ls,  k(t)  =  L(t)/Ls  (7) 


These  parameters  combine  as 

K  =  k2/[(l+k1)(l+k0)-k2](l+k0)  (8) 

From  the  ideal,  static  circuit  (R  =  dL/dt  =  0),  we 
have  the  equivalent  inductance 

LE  =  LP  [  k/(l  +  A0)  ]2/k  (9) 

The  normalized  primary  and  secondary  currents 
are,  respectively 

h  =  J,  /[Vo/(Le/C)1/2]  ,  j2  =  h  /[k1/2Vo/(Ls/C)1/2]  (10) 

The  dimensionless  parameters  for  system 

behavior  now  become: 

Dynamic  parameter 

T  =  K  [CV02/2m]/[Ls  /Ldo(LEC)1/2]2  (11a) 

Inductance  parameter 

A  =  Ls/Ldor0  (lib) 

Liner  and  circuit  resistance  parameters 

0  =  (lec)1/2/(ls/r0)  ,  ec  =  (LEQ1/2/(Ls/Rc)  (1  ic,d) 

Liner  heating  parameter 

K  =  k (prio/A2)  CV02  (LeC)1/2/Ls  (lie) 

Figures  5  displays  the  nonnalized  solutions  for 
circuit  and  liner  behavior,  with  k  =  0.9,  and 
Xo=Xi=0.1,  and  for  the  same  values  of 
parameters  used  in  the  simple  circuit. 
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Figure  5.  Normalized  solutions  for  transformer-coupled 
circuit.  Resistance-increase  q  divided  by  3.5  to  fit  graph. 

For  r0=2  cm  and  h=lcm,  and  a  12-tum  primary, 
the  initial  energy  is  460  kJ  (C=32.5  |df,  V0=167 
kV),  Ls=4  nH  and  LP=580  nH,  allowing  Li=58 
nH,  which  is  a  substantial  gain  over  the 
inductance  of  the  simple  circuit,  L0=4  nH. 
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